Finite Symmetry of Leptonic Mass Matrices 

C.S. Lam 

Department of Physics, McGill University 
Montreal, Q.C., Canada H3A 2T8 
and 

Department of Physics and Astronomy, 
University of British Columbia, 
Vancouver, BC, Canada V6T 1Z1 
Email: Lam@physics.mcgill.ca 

Abstract 

We search for possible symmetries present in the leptonic mixing data from SU(3) subgroups of 
order up to 511. Theoretical results based on symmetry are compared with global fits of experi- 



mental data in a chi-squared analysis, yielding the following results. There is no longer a group that 
can produce all the mixing data without a free parameter, but a number of them can accommodate 
the first or the second column of the mixing matrix. The only group that fits the third column 
is A(150). It predicts sin 2 2^13 = 0.11 and sin 2 2^23 = 0.94, in good agreement with experimental 
results. 
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I. INTRODUCTION 



In the days when the reactor angle of neutrino mixing was thought to be zero and the 
atmospheric angle maximal, mixings could be taken to be tri-bimaximaL and explained by 
a 1S4 symmetry without any free parameter, or groups containing it l|. Now that both 
global fits and direct measurements show the reactor angle to be non-negligible and the 



atmospheric angle possib 
explain the new data {5, 



non- maximal (2 — 4 1 , many suggestions have been advanced to 



6[. In this paper we investigate whether a finite symmetry still 
exists to accommodate them. 

The group theory of mixing reviewed in Sec. 2 will be used to extract columns of pos- 
sible mixing matrices allowed by a finite group Q. With the help of the powerful group 
software GAP [7], we examine all finite subgroups of SU(3) up to order 511 and compare 
their predictions with global fits of experimental measurements. Since phases are unknown, 
only absolute values of the mixing matrix elements are used. The experimental data and 
the global analysis used in a chi-squared comparison will be discussed in Sec. 3. The result 
presented in Sec. 4 can be summarized as follows. No group can simultaneously accommo- 
date all three columns of the mixing matrix, like S4 was able to do for the tri-bimaximal 
mixing. Many groups can accommodate the first or the second column, but the success 
is not necessarily trustworthy because those globally fitted matrix elements depend on the 
unknown CP phase. For the third column, only the group A(150) works and it predicts 
sin 26 13 = 0.11 and sin 2 2# 23 = 0.94, in good agreement with experimental data. In Sec. 5, 
after a review of the mass matrices in the presence of A4 and £4, we will present the form 
of the mass matrices for A(150). 

We close this section with a brief remark comparing the group theoretical method used 
here and other approaches. Like texture zeros, both are theories of fermion mass matrices 
within the Standard Model. There is no need to introduce additional Higgs or valons, 
together with their alignments. Clebsch-Gordan coefficients are not needed. Unlike the 
texture-zero approach, mixing parameters are determined by symmetry here and not by 
ratios of fermionic masses. The group-theory approach is consistent with dynamical models 
based on horizontal symmetry if the valon alignments are given by the invariant eigenvectors 
of the residual symmetries discussed in the next section. Such alignments can be obtained 
from group-invariant potentials in the weak-coupling approximation [8]. 
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II. GROUP THEORY OF MIXING 



Every mixing produces a Z2 x Z2 symmetry in the neutrino Majorana mass matrix and a 



Z n symmetry in the 
be reviewed below il 



eft-handed charged-lepton mass matrix. The group theory of mixing to 
"si, 9] is based on the simple assumption that these natural symmetries 
are the residual symmetries of a horizontal symmetry group. 

Let M v be the symmetric mass matrix of active neutrinos, and M e := M\M e be the 
hermitian mass matrix of left-handed charged leptons. In the basis where M e is diagonal, 
the neutrino mixing matrix U renders U T M U U diagonal. Let u\, U2, u 3 be the three columns 
of U. Then the unitary matrices defined by 

Gi = +u\u\ — U2U2 — u 3 u\, 
G 2 = —Uiu[ + u 2 u\ — u 3 ul, 

G 3 = -uiu\ - u 2 u\ + u 3 u\ [Gi] (1) 

mutually commute and commute with U. They satisfy Gf = 1, GiGj = Gk k different), 
and GjM u Gi = M v . Thus they generate aZ 2 xZ 2 symmetry of M v . In the mean time, since 
M e is diagonal, every unitary diagonal matrix F commutes with it, giving F^M e F = M e . 
If F n = 1, then F generates a symmetry group Z Tl of M e . We will assume the eigenvalues 
of F to be non-degenerate so that F diagonal forces M e to be diagonal. 

Since these symmetries are always present, additional input is required to nail down the 
mixing. To this end we shall assume F and G to be residual symmetries of some finite 
group Q, in the sense that both are members of the group. In that case the group structure 
imposes a correlation between F and G, allowing G and hence U to be determined when F 
is given in its diagonal form. We call Q a 'partial- symmetry group if it contains F and one 
Gi, and a full-symmetry group if it contains F and two mutually commuting GVs. It does 
not matter which two to choose because the third one, being the product of the first two, 
must also be in Q. 

Conversely, given a finite group Q, any of its order-2 elements is a candidate of Gi, and 
any element with order larger than 2 is a candidate of F. In a 3-dimensional irreducible 
unitary representation of Q, the invariant eigenvector (the one with eigenvalue +1) of Gi 
in the F-diagonal basis constitutes a column Ui of a possible mixing matrix U. In order to 
ensure that Gi has one +1 eigenvalue and two —1 eigenvalues as in (CQ), we shall confine to 
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those representations of Gi with character —1. Going through all combinations of F, G{, 
and 3- dimensional irreducible representations, we get all possible mixing-matrix columns 
allowed by this group. We shall refer to these theoretically predicted mixing columns as 
mixing vectors, to distinguish them from the globally fitted experimental mixing columns. 
As long as Q is a finite group, the number of allowed mixing vectors is finite, though this 
number could be large for a large group. Two such vectors U{ and Uj may fit into the same 
mixing matrix U if and only if Gi and Gj commute. 

Within this scheme, the order of Q must be even because it must contain at least one 
order-2 member Gi. If it is to be a full-symmetry group, then its order must be divisible 
by 4 because there must be two distinct order-2 elements Gi present. This last condition 
is necessary but not sufficient for a full-symmetry group also requires those three order-2 
members to mutually commute. 

Since phases of neutrino mixing have not been measured, only the absolute values of these 
columns are needed for experimental comparison. Moreover, group theory can never know 
how to label the neutrino flavor states nor the mass eigenstates, hence any mixing vector 
from group theory can be used to compare with any of the three columns of the experimental 
mixing matrix, and the rows may be permuted any way we want before the comparison. 

There is a large number of allowed mixing vectors for a large group Q to make the 
comparison a daunting task to do by hand. Fortunately, a powerful free software GAP 
is available to help us. Note however that the irreducible representations given in GAP 
may not be unitary so to use it we must first obtain from it the corresponding unitary 
representations. 

In Sec 4, we will use the recipe outlined here to obtain the allowed mixing vectors of all 
finite subgroups of SU(3) with an order less than 512. 
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III. EXPERIMENTAL DATA AND GLOBAL FITS 



14 



Recent experiments give the following values for the reactor angle #13 |4|: 

sin 2 26 13 = 0.089 ± 0.010 ± 0.005 (Daya Bay) 
sin 2 29 13 = 0.109 ± 0.030 ± 0.025 (Double Chooz) 
sin 2 2^ 13 = 0.113 ±0.013 ±0.019 (RENO) 
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sin 2 2^13 = 0.104 ±0.060 ±0.045 (T2K, normal hierarchy) 



sin 2 29 13 = 0.128 ±0.070 ±0.055 (T2K, inverted hierarchy), [reactor] 



(2) 



A preliminary result from MINOS 4J also shows that the atmospheric mixing may not be 
maximal: 



sin 2 26 23 = 0.96 ± 0.04 (MINOS, v) 

sin 2 29 23 = 0.97 ±0.03/0.08 (MINOS, V). [atm] (3) 

These results are to some extent anticipated by global fits of the data {2, 3|. The absolute 
values of their mixing-matrix elements are shown below: 
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where N and I stand for normal and inverted hierarchies, respectively, and a, b are respec- 
tively the results taken from Q and 3. Since phases are not yet measured, only absolute 
values of the matrix elements are listed and compared. Note that the 22, 23, 32, 33 matrix 
elements depend on the unknown CP phase 5, resulting in relatively large errors and may 
therefore be somewhat unreliable. 
We use the chi-square measure 

|^|) 2 /24 2 [chi2] (8) 



X 



£(l* 

•i=i 



to gauge the goodness of a theoretically predicted mixing vector c = (ci,C2,C3) T , where 
\Uij \ ± <7y is taken from one of the four globally fitted mixing matrices in ([7]). A factor of 2 
is included in the definition to simulate the two degrees of freedom in a normalized column, 
but since it is the relative size of x 2 that will be invoked, it does not matter whether we drop 
that factor or not. The result of these fits for finite subgroups of SU(3) will be discussed in 
the next section. GAP is used to produce these results, but as remarked in the last section, 
the irreducible representations given by GAP are not necessarily unitary, so they have to be 
rendered unitary first before the mixing vectors c can be computed. 

To have a standard for comparison, we list in Table I the chi-square of each column 
of the tri-bimaximal matrix. The absolute values of its third (bimaximal) column is 
(0, .707, .707) T ~ (0, 1, 1) T , that of its second (trimaximal) columns is (.577, .577, .577) T ~ 
(1, 1, 1) T , and that of its first column is (.816, .408, .408) T ~ (2, 1, 1) T . 
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Table I. The x 2 -values of the columns of a tri-bimaximal matrix 

With the newly measured reactor angle, the third column having a % 2 -value over 100 
is clearly unacceptable. The fit to the first and second columns are much more tolerable, 
but to some extent that may be due to the large errors associated with the unmeasured CP 
phase appearing in these two columns. In what follows we will reject all fits with \ 2 > 7; 
those that survive in Table I are underlined for easy comparison. 



IV. FINITE SUBGROUPS OF SU(3) 

Finite subgroups of SU(3) with a three-dimensional irreducible representation and an 



order less than 512 are tabulated 



in [l0| 



and reproduced here in Columns A and B of Table 
II. Column A gives the designation of a group in the Small Group Library of GAP; the first 
of the pair is the order of the group, and the second is the GAP-assigned number among 
groups of that order. Column B gives the popular name of the group, if there is one. If 
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the group is known under different names, then several of these may be given. Column C 
indicates whether the group contains A4 or 5*4 as a subgroup. A symbol o indicates that it 
contains A±, and a symbol • indicates that it contains S4, which then must also contain A4. 

If v4 4 is a subgroup, the group must contain the (unnormalized) trimaximal mixing vector 
(1, 1, 1) T . If S4 is a subgroup, then it must contain both (1, 1, 1) T and (2, 1, 1) T , with the 
corresponding x 2 given in Table I. Since these x 2 mav be considered as reasonable, we can 
use the group as a partial-symmetry group to build a neutrino mass matrix, with (2, 1, 1) T 
in the first column of its mixing matrix, or (1, 1, 1) T in the second column. This strategy has 
been used, for example, in refs. {(3, sj]. The corresponding mass matrices are also discussed in 
Sec. V below. However, we must not use both of them simultaneously, for if we do so then 
the group becomes a full-symmetry group and the mixing matrix is automatically TBM, 
giving rise to a zero reactor angle. 

Two other symbols also appear in Column C. The symbol x is used to indicate groups of 
odd order, which contains no element of order 2, and therefore will be ignored from now on. 
The symbol p is used to identify groups that can only be partial-symmetry groups. These 
groups do not contain two mutually commuting order-2 elements so they can never serve as 
a full-symmetry group. 

For each group of even order, we compute all its mixing vectors c = (ci, C2, c^) T using the 
recipe discussed in Sec. 2, then its x 2 ■ We reject cases where x 2 > 7 f° r an f° ur global fits. 
Otherwise the values |ci|, \cz\, [ C3 1 are listed in Column D, with the minimal x 2 among the 
four global fits appearing in Column E, and the corresponding global fit in Column F. The 
column that it fits, namely, j of |Z7y| in (jBJ), appears in column G. For groups containing A4 
as a proper subgroup, the (1, 1, 1) T mixing vector is understood and will not be listed. For 
groups containing S4 as a proper subgroup, neither the (2, 1, 1) T nor the (1, 1, 1) T appears 
explicitly. The symbol — is used to indicate that there is no fit whatsoever with x 2 < 7. 

The results appearing in Tables Ha, lib, lie can be summarized as follows. With so many 
groups and so many possible mixing vectors for each group, it is somewhat surprising that so 
few passes the experimental test. Besides the (2, 1, 1) T mixing of the first column for groups 
containing S4, and the (1, 1, 1) T mixing of the second column for groups containing A4, there 
are only a few that fit the first or the second column, albeit with a better \ 2 ■ The only 
group that really fits the third column is A(150), with a mixing vector (.170, .607, .777) T , 
which gives rise to sin 2 26 l 13 = 0.11 and sin 2 2^23 = 0.94, in good agreement with direct 
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measurements {4]. The only other group that fits the third column with a % 2 < 7 is the group 
A(294), but it yields too small a reactor angle with sin 2 2#i3 = 0.06 and sin 2 2023 — 0.97 
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Table Ha. Comparison of predictions of SU (3) subgroups with experimental data 
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Table lib. Comparison of predictions of SU (3) subgroups with experimental data 
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Table lie. Comparison of predictions of SU (3) subgroups with experimental data 

V. MASS MATRICES 

If F and G are the residual symmetries, then the invariant conditions F^M e F = M e and 
G M U G = M v can be used to determine the form of these two left-handed mass matrices. 
In this section we will illustrate this procedure by determining the mass matrices of (groups 
containing) S4, and the group A (150). 
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For S4 in the F-diagonal representation, M e = diag(a, 5, 0) is parametrized by three real 
numbers a, 5, (ft which can be used to fit the squared masses of charged leptons. It is well 
known that M v has a 2-3 and a magic symmetry parametrized by three complex parameters 
that can be used to fit the neutrino masses and Majorana phases. 

For A (150), it turns out that G is very complicated in the F-diagonal representation, 
but fairly simple in the representation given by GAP, so the mass matrices will be given in 
that representation. As a warm up, let us see what the mass matrices of S4 look like in the 
GAP-representation when 64 is taken both as a partial-symmetry group and a full-symmetry 
group. 

In the GAP representation of S4, F and Gi can be taken to be 
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Since M e is hermitian and M v is symmetric, they can be parametrized in the following way: 
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where a, 5,<p are real and the rest of the parameters are generally complex. Putting in the 
invariant conditions, it is easy to see that 7 = (3*, and (ft = 5 = a. For M e , we are therefore 
left with one real parameter a and one complex parameter j3, which can be used to fit 
the three real mass squares m^m^m^. For M u , the invariant conditions for G = G\ give 
d = a and e = — c, for G 2 give c = e = 0, and for G 3 give d = a and e = c. In each case 
when S4 is treated as a partial-symmetry group, there are four parameters left which can be 
used to determine the three neutrino masses and Majorana phases, as well as the remaining 
mixing parameters. If S4 is taken to be the full-symmetry group, when the equalities for 
Gi, G2, G3 are simultaneously satisfied, then the mixing is automatically tri-bimaximal, and 
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the neutrino mass matrix 

fa b 0\ 

M S V 4 = b a (12) 
\0 // 

is parametrized by three parameters which can be used to determine the neutrino masses 
and Majorana phases. 

For the partial-symmetry group A(150) which gives rise to the desired reactor and atmo- 
spheric angles, the residual symmetries in the GAP representation can be taken to be 
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where 5 = exp(27ii/5) is the fifth root of unity. With the invariant conditions, we can derive 
the mass matrices in the GAP representation to be 
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The charged-lepton mass matrix is identical that of S4, and the neutrino mass matrix would 
become the neutrino mass matrix of the full-mixing S4 if we put e — > and 05 — > 1. It has 
one more complex parameter e compared to the S4 case because in this case the solar angle 
and the CP phase have to be fitted. 



VI. SUMMARY 



We have used the group theory of mixing to determine whether any of the finite subgroups 
of SU(3) up to order 511 can be a symmetry group of neutrino mixing. We conclude that 
none could be a full-symmetry group, but several may serve as a partial-symmetry group 
for column 1 or column 2 of the mixing matrix. The only group where the third column 
can be accommodated is A(150), which yields sin 2 2# 13 = 0.11 and sin 2 2# 2 3 = 0.94, in good 
agreement with direct experimental measurements and global fits. The three-dimensional 
representation of A (150) where this occurs is complex, so it may be difficult to produce a 
dynamical model in this case. 
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